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FOREWORD 

This  report  ms  prepared  for  the  Applied  Mathematics  Research 
Laboratory*  Aerospace  Research  Laboratories,  by  P.R.  Krlshnalah  and 
A.K.  Chattopadhyay  under  project  7071,  'Research  In  Applied  Mathematics*. 
The  work  of  A.K.  Chattopadhyay  Is  performed  at  the  Aerospace  Research 
Laboratories  while  In  the  capacity  of  Technology  Incorporated  Visiting 
Research  Associate  under  contract  F33615-71<^-1463,  T.I.  Project  No. 
4262B. 

In  this  report,  the  authors  consider  the  problems  connected  with 
certain  non  null  distributions  associated  with  the  eigenvalues 

of  a  class  of  random  matrices. 

The  authors  wish  to  thank  Dr.  V.B.  Walkar  for  some  helpful  discus¬ 
sions.  Thanks  are  also  due  to  Mrs.  Georgene  Graves  for  typing  the  manu¬ 
script  carefully. 
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ABSTRACT 

In  this  paper,  the  authors  derived  expressions  for  the  marginal 
distributions  of  any  few  consecutive  ordered  roots,  moments  of  the 
elementary  symmetric  functions  of  the  ordered  roots  and  the  Laplace 
transformations  of  the  traces  of  a  class  of  randan  matrices  In  the 
noncentral  cases.  This  class  of  random  matrices  Includes  the  MANOYA, 
Canonical  Correlation,  and  Wlshart  matrices.  The  expressions  obtained 
here  are  In  terms  of  the  linear  combinations  of  the  products  of  double 
Integrals;  these  double  Integrals  can  be  evaluated  without  difficulty 
for  the  cases  of  the  random  matrices  that  occur  commonly  in  multivariate 
statistical  analysis.  In  deriving  the  results  presented  In  this  paper, 
the  authors  exploited  the  method  of  de  Bruljn  (J.  Indian  Hath  Soc.  19 
133-152)  for  the  evaluation  of  certain  Integral. 
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1.  INTRODUCTION 

Krishna lah  and  his  associates  [4-6]  derived  exact  expressions  for 
the  marginal  distributions  of  any  single  or  few  ordered  roots  of  a  class 
of  random  matrices  as  well  as  the  distributions  of  the  traces  of  two 
random  matrices.  They  derived  the  above  results  by  exploiting  the 
method  of  integration  over  alternate  variables  (e.g.,  see  [7]).  For 
a  very  brief  summary  of  the  literature  on  the  distributions  of  the 
Individual  roots  and  traces  of  some  random  matrices  the  reader  Is  re¬ 
ferred  to  [4-6]. 

In  general  one  can  use  any  suitable  function (s) (not  necessarily  sym¬ 
metric)  of  the  roots  of  the  appropriate  random  matrices  to  test  vari¬ 
ous  hypotheses  that  arise  on  such  problems  as  MANOVA,  canonical  cor¬ 
relation,  tests  for  equality  of  covariance  matrices,  etc.  Of  course, 
the  choice  of  these  functions  depends  upon  such  factors  as  the  optlmtm 
properties  of  the  tests  and  the  feasibility  of  the  evaluation  of  the 
distributions.  It  will  be  of  Interest  to  examine  the  distributions  of 
some  of  these  functions  which  have  at  least  Intuitive  appeal  In  testing 
some  of  the  hypotheses. 

In  this  note,  we  extend  the  results  In  [4-6]  to  the  non-central  cases 
by  using  the  method  of  de  Bruljn  [1]  for  the  evaluation  of  certain  Inte¬ 
gral.  Taking  advantage  of  de  Bruljn's  method  [1],  we  have  also  derived 
the  moments  of  the  elementary  symmetric  functions  of  roots  of  a  class  of 
random  matrices  In  the  non-central  case.  Here  we  note  that  Filial  and 
his  associates  (see  [2,3,9]  and  the  references  there)  derived  the  moments 
of  elementary  symmetric  functions  of  roots  explicitly  In  some  special  carrs. 
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2.  PRELIMINARIES 


* 

de  Bruijn  [1]  proved  the  following  uesful  result: 

Lewna  2.1  (de  Bruljn).  Let  (a,b)  be  any  Interval  finite  or  Infinite. 

Then 

|  •••  dXj ***dxn  ■  Pf(A)  (2.1) 

4  1  xi  £  *•*  £  xn  <.  b 

where 


*(x,.-.xn)  ■  [y1Jl,  yfJ  ■  *,(Xj). 

A  *  (S(j)  and  Pf(AJ  denotes  the  Pfaffian  of  A.  Here  the  elements  #1j  8re 
given  by 


a1j  " 


^(x)  aj(y)  Sgn(y  -  x)  dx  dy 


(2.2) 


1»j  ■  l»,,#»2m 

if  n  *  an;  if  n  •  an  +  1,  then  they  are  given  by  (2.2)  and  by 


a2nH-2,  2itH-2  "  °* 


b 


a1,  2m+2  *  "a2m*'2f  1 


*^(x)  dx  1  ■  !,•••, 2m+l  . 


'a 

We  need  the  following  In  the  sequel: 

Lewna  2.2.  Let  n{xj,»«*,xn)  be  a  symmetric  function  of  x^, 
♦(x-|.»*«,xn)  be  as  defined  before.  Then 


•  »  xn  and 


*  The  integral  In  Eq.  (2.2)  of  [6]  Is  a  special  case  of  the  integral  In 
Eq.  (2.1)  of  this  oaper. 
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j  •**  I  «(xr...,xn)  ♦(x1,...*xn)  dx1...dxf 

1  £  Xj  £  <  3 


*  £  -  £  n(«, )n]  {♦,<«,)  dx()  (2.3) 


where 


In  addition,  If 


"(*1 . V  ■  J)  C,  •  (2. 

*<wre  CK  depends  upon  «,.•••. k„  end  J,  denotes  the  surma t ton  over  the 
different  ..events  of  <  .  subject  to  suitable  restrictions. 


(2.3a) 


we  have 


*  i  xi  i  •••  <*n  <.  b 

•  l}  Cc  Pf(Ar) 


(2.4) 


where  Ac  *  (a^)  # 


In  the  above  lenroa,  a*j  are  given  by 
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b 


b 


(2.5) 


•M- 


K,  Kj 

x  ♦1(x)  y  J  4j(y)  Sgn(y  -  x)  dx  dy 

1,j  ■  !,•••, 2m 


when  n  »  2m;  if  n  *  2m  ♦  1,  the  elements  a*j  are  given  by  (2.5)  and 
aai*2,  2m+2  *  0 

^  an.2  ■  I  * 1  *i(x)  1  ■  ♦  '• 

Ja  i 

Lemma  2.2  can  be  proved  by  following  the  same  lines  a*  In  de  Bruijn  [1], 


3.  MOMENTS  OF  ELEMENTARY  SYMMETRIC  FUNCTIONS  OF  ROOTS 
Let  ip***»tp  be  the  latent  roots  of  a  random  matrix  and  let  us 
consider  the  situations  when  their  joint  density  In  the  non-central 
case  Is  of  the  form 

fUtt-'.iJ  ■  C  n  n  (1.  -  t.)  U  I,  a(«)  n  (L)  (3.1) 

1  p  1-1  1>J  J  2  1  * 

a  <  •••  1  *p  i  b 

where  L  *  dlag  U^»***»tp).  *  *  (kj,**«,kp)  Is  a  partition  of  k  subject 
to  suitable  restrictions,  n^L)  is  a  symmetric  function  and  a(ic)  depends 
upon  the  population  parameters  and  the  elements  of  the  partition  *. 

Also  denotes  the  sunraation  over  all  partitions  of  k,  whereas  ^ 
denotes  the  sunroatlon  over  k.  Now,  let 
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(3.2) 


r.  r. 


*  ^3  br 

m 

where  l 3  denotes  the  animation  over  r1#***frp  subject  to  some  suitable 
restrictions  and  b*  depends  on  rj,***»rp  and  k.  If  n^L)  Is  the  zonal 
polynomial  and  the  elements  of  the  partition  r  ■  (k^ » *• **kp)  are  sub¬ 
ject  to  the  restriction  k,  >_  •••  >.  k  >  0,  then  n  (L)  Is  of  the  form 

*  r  * 

(3.2). 

The  elenentary  symmetric  function  of  order  q  Is  given  by 


where 


c(t-j »*•  t1  —t^ 


i.<««  «<i 

i  q 


i4  *t)— 11( 


IctS'j  )  )S  ■  £5  f51»”**sp*)  * 


% 

’V’ 


(3.3) 


P*  ■  (J>  I5 

denotes  the  sunnatlon  over  sj,*«».sp*  such  that  s^  ♦  •••  f  sp*  * 
s,n<j  ♦  •••  ♦  np  *  qs,  and  each  n1  depends  upon  s^.'^.s^.  We  know  that 


n  n  (L)  UU,,  — .tJ}s 
?»1  1  K  1  P 


Is  a  symmetric  function  of  ip‘**tip.  So  by  (2.3)  and  (2.4)  we  get 
the  sth  moment  of  as 

,s  /b 


[UUjt^’ttp)}8}  ■  C  ^  aU) 


n  (eUj)  tj'1}  nK(L)  UUj.^'.tp})5  dV**d£p 


E(lj.*»*«tp) 


C  ^2  ^1  ^3  ^5  br(sl 


(A(r1#«*,r  ;  ^.••••np)) 


(3.4) 


where 


Mr-i J  n1.***.np)  * 


p  1+r.  +  n.,-1 

n  {*(x.)  x  dx.) 

1-1  1 


E(xlt***,xp) 


Pf(aij(r1,»«*,rp;  n1,*-*,np)) 


and 


*1J<rT— ’V  nl’"-,V 


b  ,b  iHri  *  "I-1  *  v1 


*(x)  *(y)  sgn(y  -  x)  dx  dy 

when  p  -  2m;  when  p  -  2m+lf  the  elements  a^(ri»***»rpl  V***’V  ar® 
as  given  above  and 
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leiiiwrir  Vi  ‘eTTn 


ate  tftamoiiTiiia 


*2*2,  2«H-2(rl*,**»rpl  V‘*V  “  ° 

#1,  2»*2(rr***»V  V*»V 

*  "*2m*2,  1^r1»***»rpi  ^ *np) 

P|+n>1«l 

•  J  x  *(x)  dx,  1  ■  V.2m+1. 

Th«  swnatlons  J,,  \v  \y  £g,  are  defined  In  (3.1),  (3.2)  and  (3.3). 

In  the  central  case,  let  us  assume  that  the  joint  density  of 
*•••,!  Is  of  the  form 

,  P 

D)  * c  1  *(li}  n  U«  -  O. 

»  p  i-i  1  i>j  1  j 

In  this  case 

~^UU1«#**,ip))s  •  C  ^$r***»sp*j  Pf(A{0,»«‘,0;  n1,*«*,np) 
where  A(0,***,0in^,***,np)  Is  given  by  (3.4). 

4.  DISTRIBUTION  OF  A  FEW  ORDERED  ROOTS 
We  need  the  following  In  the  sequel. 

*g£i*  A-1*  Let  (a,b)  be  any  Interval  finite  or  infinite.  Also,  let 
I(hj,‘««,kn;  a,b)  *  j  •••  |  •  (x^,*»»^n;  k^,»*»,kn)  dx^*«»dxn 

where 
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R*:  a  «  x1  <  •••  <  x,.  <^1  <  Xp+s  <Wl  i***  ixpib 

Then 

|  *•*  j  ♦(xie*»«,xp)  n(x^,*»»,Xp)  dx1***dxf  dx^^^^dXp 

*  lk*il.  iiS*-!!.  Jk  .  l 
■  Z6 17  (-')  2  ^  t-i1  1-1 

'  * 

V(Xr+i'"‘'Vs'  “•  ECxr...,xp)  E(xr+S4>i»,,**xp)  (4.3) 

R** 

V(x1(*«*,xr;  k-|,**«,lcr)  Br+s+T***,ep^ 

n(xlt-..,xp)  dXl.-dxr  <*xp+s+1  •*»dxp;  a  <  xr+1  <  —  <  Xp+S  «  b, 

and  (kj  <  <  kp)  Is  a  subset  of  (1, •••,?),  (tj  <  •••  <  tp_p)  Is  Its 

complementary  set.  Also  (a^  <  •••  <  as)  Is  a  subset  of  (t^  <  <  tp_p) 

and  (S^s+j  <  *•*  <  Bp)  Is  the  complementary  subset  of  (a^  <  •••  <  o$) 
with  respect  to  (tj  <  •••  <  tp<iip),  denotes  this  summation  over  all 
possible  (p)  choices  of  (kj,*»*,kp)  and  ][7  stand  for  the  sunmatlon  over 
(p  ~  r)  possible  choice  of  (a1  <  <  os), 

*(*»•*••■ 'V  kV" -M  ’ 


R**:  {a  ^  x1  iXp+T*  0-1, •••,»*);  x^  <.Xj  <  b,  (j  •  r+s-H, •♦•»?)). 
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i 


i 

t 


formula  (4.3)  follows  by  repeated  use  of  Laplace  expansion  of  e(x^‘»‘,xp) 
starting  with  the  first  r  columns  and  then  with  the  first  s  columns  of  the 
remaining  determinant.  Now  again  If  n(xj,**»,xp)  Is  of  the  form  n{xj,**»,xp)  * 
d.  d 

C4  x.j***xpp  such  that  n(x1»***fxp)  Is  symmetric  In  x^*»*^tp  and  Is 
similarly  defined  as  In  (2.3a) 


R* 


Is  k  h  <-’> 


♦(X'j #**•  #xp)  n(x1,**«,xp)  dx^^dXy.  dx^^^dXp 

.UJjtlU  jk,  ♦  1., 

-i)  2  2  isr  i-i1  c« 


ut  \  ^r+l  ^r+s 


r+s'l  1 


s'  r+1  r+s 


^Sr+s+l****,8p’^m+l****,dpl  xr+s*5^ 

Where  £g  Is  given  In  (2.3a)  and  I(kj,»*o,kr,d.j,»‘»,dp;  a, fa)  Is  obtained 

d 

from  I(tj,«*»,kr;  a,b)  replacing  e^fx^)  with  xf  ♦j^(xj). 

Now  let  us  take  the  joint  density  of  the  roots  In  the  noncentral  case  as 
given  fay  (3.1).  Then  applying  formula  (l.3)we  get  the  density  of  roots 

xr+r‘*',xr+s  as 

f(xr+l»"',xr+s^ 
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a(c)  b J  I(k1,---,kr,x1,-..,xr;  a.x^) 


^6r+s+l,*“,8p*ir+s+r*“'xp;  xr+s,b^ 

v(5W*’»W  ar“-*0s) 

dr+l  dr+s  . 

xr+i  ***xr+s  •  *  i  xr+l  ^  ±  xr+s  ±  b 

where 

♦<v— .*P)  ■ del  (4-5> 

C,  Ii  lz  I3  Ifi  I7  are  a^  defined  earlier  b*  a.id  are 

» 

defined  In  the  same  way  as  br,  rj,-«»,r  In  (3.2). 

5.  PROBABILITY  INTEGRALS  OF  ORDERED  ROOTS 
As  before,  lat  us  assune  that  the  joint  density  of  the  roots 
tj  <  •••  <  ip  Is  given  by  (3.1).  Then  the  density  functions  of  the 
extreme  roots  can  be  obtained  by  following  the  same  lines  as  In 
Krishnalah  and  Chang  [5]  and  using  (L',3)  and  (2.4).  Also  the  c.d.f. 
of  the  Individual  extreme  root  can  be  computed  by  applying  (2.3)  and 
(2.4).  The  probability  integral  associated  with  the  joint  density 
of  any  two  ordered  roots  Is  easily  obtained  by  following  the  same 
lines  as  In  Krishnalah  and  Walkar  [6]  and  applying  formulae  (2.3) 
and  (2.4).  The  expressions  associated  with  the  d*»ns1ty  functions 

11 


and  c.d.f.s  of  the  extreme  roots  and  the  probability  Integral  asso¬ 
ciated  with  any  two  ordered  roots  are  not  given  here  for  the  sake 
of  brevity.  We  will  now  derive  the  c.d.f.s  associated  with  the  Inter¬ 
mediate  roots. 

It  Is  known  that 

P[tr  <  x]  -  PCipj.1  <.*]  ♦  <  •••  <  ip  <  x  <  <  1^]  (5.1) 

Now  c.d.f.s  of  the  extreme  roots  can  be  evaluated  using  (2.3)  and 
(2.4).  The  second  term  In  (5.1)  can  oe  calculated  following  similar 
lines  as  In  Krlshnalah  and  Walkar  [6]  and  using  (4.3)  and  (4.4).  It 
is  given  by 


P[i.  <  ...  <  l  <  x  <  i  < 


r *  i  '  '  *r - r+1 . V 

ir(r  ♦  1}  ♦  Ik. 

C I2 1,  Z3  If  Hr  ^  •(>; 

I(k^ ,* *. ,kJ,ix^ ,». * ,xp»  a,x)  I(t^*...,tp_r,tr+1,...,xp;  x,b) 


where 


C.  Ip  ^3*  ^6*  ^1  <  ***  <  ^1  <  ***  <  *p-r^ 

and 


♦(ti.»**»tr) 

are  ss  defined  In  (4.5). 
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6.  NON  NULL  DISTRIBUTIONS  OF  TRACES 
Here  Me  use  the  results  In  [4]  and  treat  ooth  traces  discussed  there 
In  a  unified  May.  let  the  joint  density  of  the  roots  be  given  as  (3.1). 
Then 

E(e-t(»,  *  —  ♦  tp)) 

-  c  I2  I,  l3  .(<)  bj  Pf(Br) 

where 


8r  "  ^br,1j^1,j-1,***,p 


and 


>.1  j 


r.+1-l  r.+j-l  ^  *,u 

x  1  y  3  *{x)  *(y)  e-tx  e’ty  sgn(y  -  x)  dx  dy 


If  p  Is  odd,  the  Br  matrix  has  to  be  augmented  as  given  In  Section  2. 
Now  to  find  the  exact  distribution  we  have  to  Invert  this  Laplace 
transform.  To  this  end  we  note  that 


Pf(B  )  -  I  +  b  1  1  . b  i  i 

r  r  * 1 1  *  *2  •  2a-lf 


2a 


Where  the  sirmatlon  Is  over  all  possible  choices  of  l-],***,^  sub^ect 
to  this  restriction  1^  <  12,***,  <  12a  and  tbe  s^gn  +ve  or 

according  as  the  permutation  Is  even  or  odd.  Thus  Inverting  the  Laplace 
transform  we  get 


f(Y)  •  C  Si  I3  br  ^l»***,^2a**  ^ 


where  hp(1^ .••*,12a)  Is  t:.e  Inverse  Laplace  transform  of  Pf(Bp).  For 

^  m 

discussion  on  Inverse  transform  of  the  above  type  for  two  Important 
cases  we  refer  to  [4]. 
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